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Outlet Guide Vanes 


□ Engine Noise Sources 


Turbine 


Inlet Sources: 
Fan & Compressor 


Compressor 


Combustor 


Exhaust Noise: 

Fan, Jet, Turbine & Core 
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□ Rotor-Stator Source 



Acoustic Waves 



Tip Boundary Layer 


Rotor-Stator Interaction Noise 
Tonal & Broadband 
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□ Modeling Strategies 



Coupled Blade-Row Computation 




Isolated Blade-Row Computations 
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□ Some of the Technical Issues 


Coupled Blade Row Strategy (Navier-Stokes) 

> Blade/Vane Ratio Problem (Multiple-Passage Domains) 

> Information Transfer Across the Sliding Interface 

> Turbulence Modeling 

> Grid Issues (Structured v. Unstructured, Topology, Resolution) 

> Time Accuracy 

> ... 

Single Blade Row Strategy (N-S for Rotor, Euler for Stator) 

> Swirling Inflow/Outflow Type Non-Reflecting Boundary Conditions 

> Iterative Blade-Row Coupling ? 

> Grid Issues 

> Time Accuracy / Frequency Resolution 

> ... 


❖ Stringent Computational Accuracy 

> Acoustic Perturbations -0.2% of Background Flow (140 dB = 0.03 psi) 
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□ LINFLUX Tone Noise Prediction Results 

❖ Wind Tunnel Test Data 

> Realistic Configurations 

> Flow and Acoustic Data 



SDT Fan ADP Fan 1 
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Data-Theory Comparisons 
SDT Fan OGVs (3) 


Tip Speed: 
Frequency: 


7808 rpm (Approach) 
IxBPF & 2xBPF 
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SDT Cut-Off Stator (54-Vanes) SDT Cut-On Stator (26-Vanes) 


AFT Tone Power Levels: Predictions (Black), Data (Red) 


Cut-Off Stator (2xBPF) 


Cut-On Stator (IxBPF) 

Mode: (m,n) 

PWL (dB) 

Mode: (m,n) 

PWL (dB) 

(-10,0) 

113 114 

(-4,0) 

124 127 

(-10,1) 

100 100 

H,1) 

120 123 

(-10,2) 

101 106 



(-10,3) 

102 101 



Total 

114 115 

Total 

125 128 
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54-Vane Configuration: Leaned OGV (Straight) 



□ Synopsis 

• Converged TURBO and LINFLUX Solutions (Poor Quality Meanflow, “Separated” at the Hub) 

• Mixed Noise Reduction Benefits Predicted at 2xBPF (w.r.t. Radial SPLs & PWLs) 


Black: Radial OGV (Theory). 

Blue: Straight Lean OGV (Theory) 


m,n 

SPL 

PWL 

(-10,0) 

118 

111 

112 

105 

(-10,1) 

106 


100 


Total 

118 


112 
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54-Vane Configuration: Leaned OGV (Composite) 



□ Synopsis 

• Converged TURBO and LINFLUX Solutions (Meanflow Solution Could be Improved Further) 

• Sizable Noise Reduction Benefits Predicted at 2xBPF (w.r.t. Radial SPLs & PWLs) 


m,n 

SPL 

PWL 

(-10,0) 

118 

111 

112 

105 

(-10,1) 

106 


100 


Total 

118 


112 



Black: Radial OGV (Theory) 
Blue: Composite Lean (Theory) 
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Data-Theory Comparisons 
ADP Fan 1 OGV 


Tip Speed: 
Frequency: 


8750 rpm (Takeoff) 
2xBPF 
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□ Mode Power Levels 


• Highly Converged TURBO and LINFLUX Solutions 

• Excellent Data-Theory Comparisons 



Predictions (Black), Data (Red) 


Cut-Off Stator (2xBPF) 

Mode: (m,n) 

PWL 

(dB) 

(-9,0) 

122 

122 

(-9,1) 

121 

121 

(-9,2) 

119 

119 

(-9,3) 

111 

110 

Total 

126 

126 
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Upstream of OGV 


□ Mode r.m.s. Pressure Levels 
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Conclusions & Issues 


□ Need a robust mean flow solution for reliable LINFLUX results. 

□ Inviscid mean flow calculations are problematic for unconventional 
geometries. 

□ Do linearized Navier-Stokes methods offer any advantages? 

□ If so, can one do “selective” linearization? 
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Development of a High-Accuracy 
Finite-Difference, Time-Domain 
Fan Noise Prediction Code 

Ray Hixon 
R.M. Nallasamy 
Scott D. Sawyer 
Rodger W. Dyson 
Danielle Koch 
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Why High-Order Differencing? 

• In an unsteady problem, 10 , 

waves of various types 

, 7 io " 

must be propagated. 

^ 10" 1 
O 

• The errors in the numerical § io° 

W 

spatial derivatives affect the | ,,,- 
wave propagation speed. 1 


• High-order schemes allow 101 
fewer points per wavelength ,o- 
to be used. 


10 1 

Points Per Wavelength 


10 2 
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Governing Equations 

• The code is designed to solve the non-linear Euler or 
Navier-Stokes equations in 2D or 3D. 

• In Navier-Stokes mode, the code is designed to be 
either a DNS solver (no turbulence model), a LES 
solver (constant-coefficient Smagorinsky subgrid 
model), or an unsteady RANS solver (with a k-s 
turbulence model). 
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Code Structure 


• The code solves the flow equations in chain 
rule curvilinear form (non-conservative). 

• The code is written in Fortran 90 with MPI 
message passing for computational efficiency, 
and is designed to be fully portable between 
computer architectures and operating systems 
(testing is currently performed on SGI, Linux, 
and Mac OSX). 

• The code uses structured multi-block grids. 
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Solution Procedure 


• The code uses finite-differences to obtain the spatial 
derivatives (explicit 2 nd order, explicit 6 th order, 7-point 
DRP, or compact 6 th order derivatives are 
implemented). 

• The code marches explicitly in time, using an 
optimized Runge-Kutta scheme. In future, a fourth- 
order Adams-Bashforth scheme will be implemented. 

• The code currently uses constant-coefficient 10 th 
order artificial dissipation. 
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Assessment 


• In previous work with an earlier version of this code, 
the benchmark problem of the gust response of a 
Joukowski airfoil was solved. 

• This test case evaluated the ability of the code to 
capture the effects of changing the airfoil geometry, 
the gust geometry, and the gust reduced frequency. 
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Closeup of Cambered Airfoil Grid 



In this benchmark CAA problem, 
the effects of wall geometry, 
gust geometry, curvilinear grids, 
and farfield boundary conditions 
are tested. 


k = 1.0 
ID gust 


k = 1.0 
2D gust 



Perturbation Pressure 
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1-DGust k = 0.1 1-DGust. k= 1.0 2-D Gust, k = 0.1 2-D Gust, k = 1 .0 


Airfoil 






Cambered 

Airfoil 






GUST3D Results 
Computed Results 
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1-DGust k = 0.1 1-DGust. k= 1.0 2-D Gust, k = 0.1 2-D Gust, k = 1 .0 


Symmetric 

Airfoil 




Cambered 

Airfoil 






R = 1 

R = 3 

R = 2 

ii 

DC 


GUST3D Results 
Computed Results 
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Future Directions 


The code is currently being parallelized. 

New boundary conditions are being added to the 
code. 

Plan to include improved artificial dissipation 
models, time stepping methods, and 
parallelization techniques. 
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2D Cascade Benchmark Test 


Alternative High-Order Approaches 


Dr. Rodger W. Dyson 
Rodger.W.Dyson@grc.nasa.gov 
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Cascade Benchmark Problem 


• Gust - Cascade Interaction Problem 

• Periodicity Requires 27 Passages (22B / 54V) 

• Gust has a Multi-Frequencies Character 

- lx, 2x&3xBPF 

- Amplitudes - 9%, 0.9% & 0.2% of the Mean Velocity 

- Minimum Wavelength is on order of 3/1 1 of the Chord 

• Accuracy Requirement: ~1% Error at 3xBPF 


Need 6 th Order Accuracy in Space & Time 
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Team Effort 


• Dr. R. Hixon, Principal Code Designer 

• Dr. R. Nallasamy, Boundary Conditions 

• Dr. S. Sawyer, Boundary Conditions 

• Ms. D. Koch, PE, Grid Generation 

• Dr. R. Dyson, Team Coordinator 

• Dr. E. Envia, Turbomachinery 
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Preliminary Cascade Results 

• The grid used by the code for this case has a 6-way 
grid singularity upstream of the leading edge. 

• Initial results are promising for this case. 
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Arbitrary High-Order Methods 

• Motivation : 

- High Resolution and Efficiency 

• Challenges : 

- Need High Accuracy in Space and Time 

- Consistent Boundary Conditions (Surface & Farfield) 

- Complex Geometry (Cartesian vs. Curvilinear) 
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Wjiven^mh^r Error Magnitude 


Why Arbitrarily High-Order? 



I st Order 
3rd Order 
5th Order 
7th Order 
9th Order 

I I th Order 
1 3th Order 
1 5 th Order 
1 7th Order 
1 9 th Order 
21 st Order 



I st Order 
— — — - 3rd Order 

5th Order 
7th Order 
9th Order 

I I th Order 
1 3th Order 
1 5th Order 
17th Order 
1 9th Order 
21 st O rder 


5 10 15 20 

-Log 10 {Error Level Required) 
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Consistent Boundary Conditions 


Propagating waves accurately in time: 


dt 


dr 


2 ! 


Requires high order time derivatives 
Otherwise will get dispersion/dissipation 

r 


dp 

~dt 


u 




dp 

dx 


dp \ du 

v 1- yp 


dy 


dx 


W 




dx 


dyj) 
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Complex Geometry - Cartesian Grid 

• Advantages 

■ No metrics 

■ No singularities 

■ Easy grid generation 

■ Efficiency (few boundary pts) 

• Challenges 

■ Surface interpolation algorithm 

■ Resolving curvature 

■ Adaptive resolution with h and p refinement 
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Complex Geometry - Curvilinear Grid 



- 0.5 

X 


13.86 


- 0.475 


13.84 


Advantages 

■ Easy interpolation 

■ Curvature more easily resolved 

13.82 

■ Centered boundary stencils with ghost points 


13.8 

• Challenges 

- 0.525 

■ Computing very high-order metric derivatives 

■ 1 st order grid singularities 

■ High order boundary conditions are more complex 


1 S2.728 

143 213 

1 23.699 

104.184 

84.6600 

65.1554 

45.641 

26.1265 

6.61204 

12.9024 

32.4169 

51.9314 

71.445R 
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Future Work 


• Validate Current compact 6 th Order Code 

- 2D benchmark problem 

• Incorporate New Technology as Needed 

- High order boundary conditions 

- Higher order time advancement everywhere 

• Validate Full 3D-Stator 

- Assess the overall efficiency/usefulness 
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Computational AeroAcoustics for 
Fan Noise Prediction 


Ed Envia Ray Hixon Rodger Dyson 

NASA GRC University of Toledo NASA GRC 


An overview of the current state-of-the-art in computational aeroacoustics as applied to fan 
noise prediction at NASA Glenn is presented. Results from recent modeling efforts using three- 
dimensional inviscid formulations in both frequency and time domains are summarized. In 
particular, the application of a frequency-domain method, called LINFLUX, to the computation of 
rotor-stator interaction tone noise is reviewed and the influence of the background inviscid flow on 
the acoustic results is analyzed. It has been shown that the noise levels are very sensitive to the 
gradients of the mean flow near the surface and that the correct computation of these gradients for 
highly loaded airfoils is especially problematic using an inviscid formulation. The ongoing 
development of a finite-difference time-marching code that is based on a b^-order compact scheme 
is also reviewed. Preliminary results from the nonlinear computation of a gust-airfoil interaction 
model problem demonstrate the fidelity and accuracy of this approach. Spatial and temporal features 
of the code as well as its multi-block nature are discussed. Finally, latest results from an ongoing 
effort in the area of arbitrarily high-order methods are reviewed and technical challenges associated 
with implementing correct high-order boundary conditions are discussed and possible strategies for 
addressing these challenges are outlined. 



